Observation-Corrected Precipitation Estimates in GEOS-5

Several GEOS-5 applications, including the GEOS-5 seasonal forecasting system and the MERRA-Land data product, rely on global precipitation data that have been corrected with satellite and or gauge-based precipitation observations. This document describes the methodology used to generate the corrected precipitation estimates and their use in GEOS-5 applications. The corrected precipitation estimates are derived by disaggregating publicly available, observationally based, global precipitation products from daily or pentad totals to hourly accumulations using background precipitation estimates from the GEOS-5 atmospheric data assimilation system. Depending on the specific combination of the observational precipitation product and the GEOS-5 background estimates, the observational product may also be downscaled in space. The resulting corrected precipitation data product is at the finer temporal and spatial resolution of the GEOS-5 background and matches the observed precipitation at the coarser scale of the observational product, separately for each day (or pentad) and each grid cell

Spring Hydrology Determines Summer Net Carbon Uptake in Northern Ecosystems

Increased photosynthetic activity and enhanced seasonal CO2 exchange of northern ecosystems have been observed from a variety of sources including satellite vegetation indices (such as the Normalized Difference Vegetation Index; NDVI) and atmospheric CO2 measurements. Most of these changes have been attributed to strong warming trends in the northern high latitudes (greater than or equal to 50N). Here we analyze the interannual variation of summer net carbon uptake derived from atmospheric CO2 measurements and satellite NDVI in relation to surface meteorology from regional observational records. We find that increases in spring precipitation and snow pack promote summer net carbon uptake of northern ecosystems independent of air temperature effects. However, satellite NDVI measurements still show an overall benefit of summer photosynthetic activity from regional warming and limited impact of spring precipitation. This discrepancy is attributed to a similar response of photosynthesis and respiration to warming and thus reduced sensitivity of net ecosystem carbon uptake to temperature. Further analysis of boreal tower eddy covariance CO2 flux measurements indicates that summer net carbon uptake is positively correlated with early growing-season surface soil moisture, which is also strongly affected by spring precipitation and snow pack based on analysis of satellite soil moisture retrievals. This is attributed to strong regulation of spring hydrology on soil respiration in relatively wet boreal and arctic ecosystems. These results document the important role of spring hydrology in determining summer net carbon uptake and contrast with prevailing assumptions of dominant cold temperature limitations to high-latitude ecosystems. Our results indicate potentially stronger coupling of boreal/arctic water and carbon cycles with continued regional warming trends

Variational assimilation of remote sensing data for land surface hydrologic applications

Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Civil and Environmental Engineering, 2000.Includes bibliographical references (p. 283-192).Soil moisture plays a major role in the global hydrologic cycle. Most importantly, soil moisture controls the partitioning of available energy at the land surface into latent and sensible heat fluxes. We investigate the feasibility of estimating large-scale soil moisture profiles and related land surface variables from low-frequency (L-band) passive microwave remote sensing observations using weak-constraint variational data assimilation. We extend the iterated indirect representer method, which is based on the adjoint of the hydrologic model, to suit our application. The four-dimensional (space and time) data assimilation algorithm takes into account model and measurement uncertainties and provides optimal estimates by implicitly propagating the full error covariances. Explicit expressions for the posterior error covariances are also derived. We achieve a dynamically consistent interpolation and extrapolation of the remote sensing data in space and time, or equivalently, a continuous update of the model predictions from the data. Our hydrologic model of water and energy exchange at the land surface is expressly designed for data assimilation. It captures the key physical processes while remaining computationally efficient. The assimilation algorithm is tested with a series of experiments using synthetically generated system and measurement noise. In a realistic environment based on the Southern Great Plains 1997 (SGP97) hydrology experiment, we assess the performance of the algorithm under ideal and non ideal assimilation conditions. Specifically, we address five topics which are crucial to the design of an operational soil moisture assimilation system. (1) We show that soil moisture can be satisfactorily estimated at scales finer than the resolution of the brightness images (downscaling), provided sufficiently accurate fine-scale model inputs are available. (2) The satellite repeat cycle should be shorter than the average interstorm period. (3) The loss of optimality by using shorter assimilation intervals is offset by a substantial gain in computational efficiency. (4) Soil moisture can be satisfactorily estimated even if quantitative precipitation data are not available. (5) The assimilation algorithm is only weakly sensitive to inaccurate specification of the soil hydraulic properties. In summary, we demonstrate the feasibility of large-scale land surface data assimilation from passive microwave observations.by Rolf H. Reichle.Ph.D

Inferring Land Surface Model Parameters for the Assimilation of Satellite-Based L-Band Brightness Temperature Observations into a Soil Moisture Analysis System

The Soil Moisture and Ocean Salinity (SMOS) satellite mission provides global measurements of L-band brightness temperatures at horizontal and vertical polarization and a variety of incidence angles that are sensitive to moisture and temperature conditions in the top few centimeters of the soil. These L-band observations can therefore be assimilated into a land surface model to obtain surface and root zone soil moisture estimates. As part of the observation operator, such an assimilation system requires a radiative transfer model (RTM) that converts geophysical fields (including soil moisture and soil temperature) into modeled L-band brightness temperatures. At the global scale, the RTM parameters and the climatological soil moisture conditions are still poorly known. Using look-up tables from the literature to estimate the RTM parameters usually results in modeled L-band brightness temperatures that are strongly biased against the SMOS observations, with biases varying regionally and seasonally. Such biases must be addressed within the land data assimilation system. In this presentation, the estimation of the RTM parameters is discussed for the NASA GEOS-5 land data assimilation system, which is based on the ensemble Kalman filter (EnKF) and the Catchment land surface model. In the GEOS-5 land data assimilation system, soil moisture and brightness temperature biases are addressed in three stages. First, the global soil properties and soil hydraulic parameters that are used in the Catchment model were revised to minimize the bias in the modeled soil moisture, as verified against available in situ soil moisture measurements. Second, key parameters of the "tau-omega" RTM were calibrated prior to data assimilation using an objective function that minimizes the climatological differences between the modeled L-band brightness temperatures and the corresponding SMOS observations. Calibrated parameters include soil roughness parameters, vegetation structure parameters, and the single scattering albedo. After this climatological calibration, the modeling system can provide L-band brightness temperatures with a global mean absolute bias of less than 10K against SMOS observations, across multiple incidence angles and for horizontal and vertical polarization. Third, seasonal and regional variations in the residual biases are addressed by estimating the vegetation optical depth through state augmentation during the assimilation of the L-band brightness temperatures. This strategy, tested here with SMOS data, is part of the baseline approach for the Level 4 Surface and Root Zone Soil Moisture data product from the planned Soil Moisture Active Passive (SMAP) satellite mission

Spring hydrology determines summer net carbon uptake in northern ecosystems

Increased photosynthetic activity and enhanced seasonal CO2 exchange of northern ecosystems have been observed from a variety of sources including satellite vegetation indices (such as the normalized difference vegetation index; NDVI) and atmospheric CO2 measurements. Most of these changes have been attributed to strong warming trends in the northern high latitudes (50° N). Here we analyze the interannual variation of summer net carbon uptake derived from atmospheric CO2 measurements and satellite NDVI in relation to surface meteorology from regional observational records. We find that increases in spring precipitation and snow pack promote summer net carbon uptake of northern ecosystems independent of air temperature effects. However, satellite NDVI measurements still show an overall benefit of summer photosynthetic activity from regional warming and limited impact of spring precipitation. This discrepancy is attributed to a similar response of photosynthesis and respiration to warming and thus reduced sensitivity of net ecosystem carbon uptake to temperature. Further analysis of boreal tower eddy covariance CO2 flux measurements indicates that summer net carbon uptake is positively correlated with early growing-season surface soil moisture, which is also strongly affected by spring precipitation and snow pack based on analysis of satellite soil moisture retrievals. This is attributed to strong regulation of spring hydrology on soil respiration in relatively wet boreal and arctic ecosystems. These results document the important role of spring hydrology in determining summer net carbon uptake and contrast with prevailing assumptions of dominant cold temperature limitations to high-latitude ecosystems. Our results indicate potentially stronger coupling of boreal/arctic water and carbon cycles with continued regional warming trends

Assessment of MERRA-2 Land Surface Energy Flux Estimates

In the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) system the land is forced by replacing the model-generated precipitation with observed precipitation before it reaches the surface. This approach is motivated by the expectation that the resultant improvements in soil moisture will lead to improved land surface latent heating (LH). Here we assess aspects of the MERRA-2 land surface energy budget and 2 m air temperatures (T(sup 2m)). For global land annual averages, MERRA-2 appears to overestimate the LH (by 5 W/sq m), the sensible heating (by 6 W/sq m), and the downwelling shortwave radiation (by 14 W/sq m), while underestimating the downwelling and upwelling (absolute) longwave radiation (by 10-15 W/sq m each). These results differ only slightly from those for NASA's previous reanalysis, MERRA. Comparison to various gridded reference data sets over Boreal summer (June-July-August) suggests that MERRA-2 has particularly large positive biases (>20 W/sq m) where LH is energy-limited, and that these biases are associated with evaporative fraction biases rather than radiation biases. For time series of monthly means during Boreal summer, the globally averaged anomaly correlations (R(sub anom)) with reference data were improved from MERRA to MERRA-2, for LH (from 0.39 to 0.48 vs. GLEAM data) and the daily maximum T(sup 2m) (from 0.69 to 0.75 vs. CRU data). In regions where T(sup 2m) is particularly sensitive to the precipitation corrections (including the central US, the Sahel, and parts of south Asia), the changes in the T(sup 2m) R(sub anom) are relatively large, suggesting that the observed precipitation influenced the T(sup 2m) performance

Uncertainty Quantification of GEOS-5 L-band Radiative Transfer Model Parameters Using Bayesian Inference and SMOS Observations

Uncertainties in L-band (1.4 GHz) radiative transfer modeling (RTM) affect the simulation of brightness temperatures (Tb) over land and the inversion of satellite-observed Tb into soil moisture retrievals. In particular, accurate estimates of the microwave soil roughness, vegetation opacity and scattering albedo for large-scale applications are difficult to obtain from field studies and often lack an uncertainty estimate. Here, a Markov Chain Monte Carlo (MCMC) simulation method is used to determine satellite-scale estimates of RTM parameters and their posterior uncertainty by minimizing the misfit between long-term averages and standard deviations of simulated and observed Tb at a range of incidence angles, at horizontal and vertical polarization, and for morning and evening overpasses. Tb simulations are generated with the Goddard Earth Observing System (GEOS-5) and confronted with Tb observations from the Soil Moisture Ocean Salinity (SMOS) mission. The MCMC algorithm suggests that the relative uncertainty of the RTM parameter estimates is typically less than 25 of the maximum a posteriori density (MAP) parameter value. Furthermore, the actual root-mean-square-differences in long-term Tb averages and standard deviations are found consistent with the respective estimated total simulation and observation error standard deviations of m3.1K and s2.4K. It is also shown that the MAP parameter values estimated through MCMC simulation are in close agreement with those obtained with Particle Swarm Optimization (PSO)

A Data-Driven Approach for Daily Real-Time Estimates and Forecasts of Near-Surface Soil Moisture

NASAs Soil Moisture Active Passive (SMAP) mission provides global surface soil moisture retrievals with a revisit time of 2-3 days and a latency of 24 hours. Here, to enhance the utility of the SMAP data, we present an approach for improving real-time soil moisture estimates (nowcasts) and for forecasting soil moisture several days into the future. The approach, which involves using an estimate of loss processes (evaporation and drainage) and precipitation to evolve the most recent SMAP retrieval forward in time, is evaluated against subsequent SMAP retrievals themselves. The nowcast accuracy over the continental United States (CONUS) is shown to be markedly higher than that achieved with the simple yet common persistence approach. The accuracy of soil moisture forecasts, which rely on precipitation forecasts rather than on precipitation measurements, is reduced relative to nowcast accuracy but is still significantly higher than that obtained through persistence

Assessment and enhancement of MERRA land surface hydrology estimates

The Modern-Era Retrospective Analysis for Research and Applications (MERRA) is a state-of-the-art reanalysis that provides, in addition to atmospheric fields, global estimates of soil moisture, latent heat flux, snow, and runoff for 1979 present. This study introduces a supplemental and improved set of land surface hydrological fields ("MERRA-Land") generated by rerunning a revised version of the land component of the MERRA system. Specifically, the MERRA-Land estimates benefit from corrections to the precipitation forcing with the Global Precipitation Climatology Project pentad product (version 2.1) and from revised parameter values in the rainfall interception model, changes that effectively correct for known limitations in the MERRA surface meteorological forcings. The skill (defined as the correlation coefficient of the anomaly time series) in land surface hydrological fields from MERRA and MERRA-Land is assessed here against observations and compared to the skill of the state-of-the-art ECMWF Re-Analysis-Interim (ERA-I). MERRA-Land and ERA-I root zone soil moisture skills (against in situ observations at 85 U.S. stations) are comparable and significantly greater than that of MERRA. Throughout the Northern Hemisphere, MERRA and MERRA-Land agree reasonably well with in situ snow depth measurements (from 583 stations) and with snow water equivalent from an independent analysis. Runoff skill (against naturalized stream flow observations from 18 U.S. basins) of MERRA and MERRA-Land is typically higher than that of ERA-I. With a few exceptions, the MERRA-Land data appear more accurate than the original MERRA estimates and are thus recommended for those interested in using MERRA output for land surface hydrological studies

Soil Moisture Data Assimilation

Accurate knowledge of soil moisture at the continental scale is important for improving predictions of weather, agricultural productivity and natural hazards, but observations of soil moisture at such scales are limited to indirect measurements, either obtained through satellite remote sensing or from meteorological networks. Land surface models simulate soil moisture processes, using observation-based meteorological forcing data, and auxiliary information about soil, terrain and vegetation characteristics. Enhanced estimates of soil moisture and other land surface variables, along with their uncertainty, can be obtained by assimilating observations of soil moisture into land surface models. These assimilation results are of direct relevance for the initialization of hydro-meteorological ensemble forecasting systems. The success of the assimilation depends on the choice of the assimilation technique, the nature of the model and the assimilated observations, and, most importantly, the characterization of model and observation error. Systematic differences between satellite-based microwave observations or satellite-retrieved soil moisture and their simulated counterparts require special attention. Other challenges include inferring root-zone soil moisture information from observations that pertain to a shallow surface soil layer, propagating information to unobserved areas and downscaling of coarse information to finer-scale soil moisture estimates. This chapter summarizes state-of-the-art solutions to these issues with conceptual data assimilation examples, using techniques ranging from simplified optimal interpolation to spatial ensemble Kalman filtering. In addition, operational soil moisture assimilation systems are discussed that support numerical weather prediction at ECMWF and provide value-added soil moisture products for the NASA Soil Moisture Active Passive mission